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Abstract 

It is shown that with allowance for nonadditivity the (canonical) 
distribution for a black hole transforms to a (microcanonical) dis- 
tribution which depends on nonadditive integrals of the motion and 
represents a nonlinear function of energy - the square of black hole 
mass. From a comparison with the standard microcanonical distri- 
bution it is found that the statistical weight of black hole should be 
proportional to the black hole area. A definition of the black hole 
entropy consistent with the statistical weight is proposed. 



The statistical interpretation of the thermodynamical entropy of a black 
hole 

Sbh = (1) 

where A is the area of black hole horizon and lp the Planck length, is still 
a central problem in black hole physics. In statistical mechanics all proper- 
ties of a system are encoded in its distribution function [1]. For a classical 
system the distribution function p(qi,Pi) determines the probability to find 
the coordinates and momenta pi of the system in intervals between 
Pi and qi + dqi, pi + dpi. For a quantum system the distribution function 
w n determines the probability to find the system in a state with energy E n . 
The determination of this function for a system is the fundamental problem 
of statistical physics. The form of the function is usually postulated; its 
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justification lies in the agreement between results derived from it and the 
thermodynamic properties of a system. 

There exist three main forms of the distribution function which play an 
important part in statistical physics, corresponding in turn to a microcanon- 
ical ensemble, whose elements have the same energy, volume and number of 
particles, a canonical ensemble, with elements of same temperature, volume 
and number of particles and a grand canonical ensemble, a set of elements 
having equal temperature, volumes and chemical potentials. The canonical 
ensemble is related to the microcanonical ensemble in that a large system 
with finite interactions can be subdivided into a relatively small independent 
subsystems; the system as a whole is isolated and is described by a micro- 
canonical ensemble. In turn, the canonical ensemble can be viewed as a "mi- 
crocanonical" with respect to grand canonical 1 . Although the concept of the 
distribution function is successfully used in statistical physics, the important 
fact is that arguments used in deriving the distribution function are based 
on the additivity properties of ordinary matter and, more fundamentally, 
on the possibility of describing a given system as made up of independent 
subsystems. However the black holes are not conventional systems: they are 
nonadditive. The main reason for this depart from the conventional systems 
lies in the long-range behavior of gravitational forces. As is well known, the 
fact that the gravitational energy is nonadditive appears already in Newto- 
nian gravity. In general relativity a local definition of mass is not possible; 
the ADM and Komar definitions of mass demonstrate this very clear. More- 
over the black hole entropy (1) goes as the square of mass M 2 in a sharp 
contrast with the additive character of entropy in ordinary thermodynamics. 
It is clear that this aspect of nonadditivity can not be ignored in deriving 
the distribution function for a black hole as well as in deriving the black hole 
entropy. But this aspect has not yet received attention in the literature. 

In this paper I suggest the forms of distribution function and entropy for 
a Schwarzschild black hole making allowance of nonadditivity. 

We begin with conventional systems. The standard determination of 
the distribution function for them is given in detail by Landau and Lifshitz 
[1]. Following Landau and Lifshitz consider an ordinary isolated system 

^dn this paper we use only the microcanonical and canonical ensembles. For this reason 
we shall not discuss the interesting approach of Gour (treating the Schwarzschild black 
hole as a grand canonical ensemble) and the reader is referred to the original paper [2]. 
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consisting of quasi-isolated subsystems in thermal equilibrium. According to 
Liouville's theorem the distribution function of isolated system is an integral 
of the motion. Due to the statistical independence of subsystems and, as a 
consequence, multiplicativity of their distribution functions, the logarithm of 
the distribution function must be not merely an integral of the motion, but an 
additive integral of the motion. There exist only seven independent additive 
integrals of the motion: the energy, the three components of the momentum 
vector and the three components of the angular momentum vector. It can be 
shown that the statistical state of a system executing a given motion depends 
only on its energy. Thus we can deduce that the logarithm of the distribution 
function must be a linear function of its energy of the form 

InwW = «W + /3E£\ (2) 

with constant coefficients a and /3, of which a is the normalization constant 
and f3 must be the same for all subsystems in a given isolated system; the 
suffix "a" refers to the subsystem a. Note that assuming another depen- 
dence lnty on i? we may not obtain an additive function on the right side 
of (2); for example, E 2 is already a nonadditive function. Since the values 
of nonadditive integrals of the motion do not affect the statistical properties 
of ordinary isolated system, these properties can be described by any func- 
tion which depends only on the values of the additive integrals of the motion 
for the system and which satisfies Liouville's theorem. The simplest such 
function is 

dw = const x 5(E - E ) JJ dT a , (3) 

a 

where the number of states of the whole system dT is a product dT = Y[ dT a 

a 

of the numbers dT a of the subsystems (such that the sum of the energies of 
the subsystems lies in the interval of energy of the whole system dE). It 
defines the probability of finding the system in any of the dT states; dT plays 
a role analogous to the classical elementary volume dq dp. The factor const 
is the normalization constant, S(E — E ) the delta function. The distribution 
(3) is called microcanonical. Note that (2) is nothing but the canonical 
distribution if we identify (3 = — l/T, a = F/T, F being the free energy and 
T the temperature of the system. 

As is easily seen, the statistical independence and additivity play a crucial 
role in deriving the distribution function for the conventional systems. Now 
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consider the black holes. An isolated black hole is not stable in empty flat 
space. So it is helpful to begin by placing black holes in an appropriate "box" 
so that a stable equilibrium configuration of black holes and their radiation 
can exist. As is well known from the thermodynamical considerations [3], at 
sufficiently high energies the entropy of such an isolated system is dominated 
by a single black hole; this corresponds to a microcanonical ensemble. As 
mentioned above, the black holes are not conventional systems: they possess 
the nonadditive properties. Because of the nonadditivity, the black holes can 
not be thought as made up of any independent subsystems. The presence of 
the logarithm in (2) is just determined by the statistical independence of the 
subsystems. Therefore, if we want to establish the distribution function for 
the black holes, we should remove the restrictions of the statistical indepen- 
dence and additivity of integrals of the motion for the subsystems of black 
hole. So dropping the logarithm and suffix "a" in (2) we obtain 

w n = f(E n ), (4) 

where f(E n ) represents a nonadditive integral of the motion and is a non- 
linear function of the black hole energy. The simplest such function is the 
square of energy, f(E) = ^E 2 . Note that as in ordinary statistics, the form of 
the distribution function must be regarded only as a postulate, to be justified 
solely on the basis of agreement of its predictions with the thermodynam- 
ical properties of black holes. Our considerations are intended to make it 
plausible, and nothing more. As a result, the (canonical) distribution for a 
subsystem (2) transforms to the (microcanonical) distribution for the whole 
system of the form 

w n = lEl (5) 

where 7 is a constant coefficient. Similarly dropping the product and suffix 
"a" in (3), we obtain 

dw = const x 5(E — E )dT, (6) 

or after integration, 

w = const x Ar, (7) 

where Ar is already the total number of states that accessible to the whole 
system in a given state. The functions (5) and (7) are obviously the same 
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and satisfy the same normalization condition J2 w n — 1- So comparing (5) 

n 

and (7) we get 

Ar oc M 2 , (8) 

where the energy of black hole is identified with its ADM mass, M. Thus 
the statistical weight of black hole should be of the form 

as is evident from dimensionality considerations. Here £ is a new constant 
coefficient, £ > 1. 

Now we can define the entropy of black hole. The entropy plays a par- 
ticularly fundamental role when the microcanonical ensemble is used. The 
entropy of a conventional system is given by the Boltzmann formula 

S = lnAr. (10) 

So assuming the usual interpretation of entropy we would have to take a 
logarithm of (9). It appears that in this case the generalized second law of 
black hole thermodynamics would be violated [4]. The argument involves 
a well-known example with the collision of black holes: two identical black 
holes collide, merge, radiating gravitational wave energy, and form a third 
black hole. According to (10), the initial entropy of the system is 

5 i = 21nAr i = 21n^A^. (n) 

On the one hand, the final entropy is bounded from above by 

S/ =lnr / = l„(i|). (12) 

On the other hand, by virtue of the generalized second law it must be greater 
then initial entropy. So we have 

21,1 <s<in S)- < i3 » 
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As is easily seen, these inequalities are satisfied only for A < 4£/f>. This 
means that the standard interpretation of the entropy in terms of the loga- 
rithm of Ar violates the generalized second law. Thus we conclude that the 
statistical interpretation of the Bekenstein-Hawking entropy is correct only 
if log is deleted from the Boltzmann formula (10), that is, the black hole 
entropy should be proportional to the number of states, 

Sbh ~ Ar. (14) 

In [4] a new statistical interpretation of the black hole entropy was suggested 
which does not use the properties of additivity and statistical independence. 
Namely, in [4] it was shown that the black hole entropy is of the form 

S BH = 27rAr, (15) 

that is, in contrast to the entropy of ordinary matter, without the logarithm, 
and the number of microstates for a given area is Ar = A/87rlp. This agrees 
with our form of the distribution function so we can identify £ = 8n in (9). 
Since our distribution is determined by the nonadditive integral of the motion 
(5), we consider the microcanonical entropy Sbh with fixed area rather than 
fixed energy. 

As is well known, a description of a black hole via a thermal ensemble 
is inappropriate. Thermodynamically the heat capacity of a Schwarzschild 
black hole is negative, so it cannot come to equilibrium with an infinite heat 
bath. From the point of view of statistical physics this could be explained 
as follows. According to the canonical distribution, the probability Pi of a 
system being in a state of energy Ei is proportional to the Boltzmann factor, 

Pi ~Ar,exp(-^), (16) 

where T is the temperature of the system. Assuming the usual interpretation 
of entropy (10), the statistical weight of black hole should grow with M as 

exp(47rGM 2 ). (17) 

In that case, the total probability diverges. However, in the case where the 
entropy is given by (15), the statistical weight grows as 

2GM 2 , (18) 
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and the probability can converge. So, it seems that the usual value of the 
statistical weight (17) better explains the breakdown of the canonical ensem- 
ble for black holes than suggested (18). But this is not the case. First, it 
is clear that an infinite heat bath is gravitationally unstable. On the other 
hand, there is always a size of bath at which the interaction energy between 
the members of ensemble is not negligible. In both cases the the canonical 
distribution is not applicable. Secondly, a black hole possesses a very special 
property which singles it out; namely, its size and temperature are not in- 
dependent parameters. As a result, the temperature of black hole does not 
remain constant, contrary to the definition of canonical ensemble. Note that 
the ordinary self-gravitating systems, for example, stars and galaxies, also 
have the negative heat capacity. And although their statistical weights grow 
not so fast as (17), they can not be in a thermal equilibrium with a heat 
bath. Thus the apparent divergence of canonical distribution, in the case 
where the statistical weight grows as (17), cannot be an evidence in support 
of the usual interpretation of black hole entropy. 
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